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AVANT-PROPOS : Le présent document constitue le chapitre 8 du traité Modélisation numérique 
discrête des matériaux granulaires, ouvrage collectif, sous la direction de Farhang Radjaï et Frédéric 
Dubois, publié dans la collection « mécanique et ingéniérie des matériaux » (MIM) aux éditions Lavoisier 
en 2010. Le titre, la pagination et l'indexation des références en sont différents. 

Il dresse un bilan de l'influence des différents paramètres micromécaniques que l'on introduit dans la mo- 
délisation des matériaux granulaires à l'échelle des grains et des contacts sur le comportement mécanique 
macroscopique en déformation quasi-statique ou en écoulement. On insiste en particulier sur la définition 
de paramètres de contrôle sans dimension, qui permettent une classification commode et générique de 
différents régimes rhéologiques. 

Ce texte contient quelques références aux autres chapitres du même traité mais peut se lire indépendam- 
ment. 

1 Introduction 

La simulation d'un matériau granulaire modélisé à l'échelle de ses constituants microsco- 
piques fait intervenir de nombreux paramètres, liés à la description géométrique du système, aux 
lois d'interaction, et aux sollicitations extérieures. 

Dans ce chapitre nous tentons un bilan de l'influence de ces paramètres sur le comporte- 
ment mécanique de matériaux granulaires en nous limitant aux modèles les plus simples, avec 
des grains sphériques ou circulaires, et en privilégiant les situations qui permettent de détermi- 
ner des comportements mécaniques exprimés par des lois constitutives. Ces lois sont liées à des 
mécanismes microscopiques que nous évoquons brièvement. Les comportements observés en si- 
mulation numérique peuvent se comparer quantitativement aux expériences de laboratoire avec 
certains matériaux modèles - assemblages de billes de verre par exemple. Avec les matériaux 
plus complexes, naturels ou manufacturés, qui intéressent directement l'ingénieur ou le physicien 
appliqué, les comparaisons sont plus qualitatives. Toutefois ces « vrais » matériaux restent régis 
par des lois analogues, issues de mécanismes rhéophysiques similaires, pour lesquelles l'étude de 
systèmes modèles fournit d'utiles indications d'ordres de grandeur. 

La partie [2j à partir des modèles micromécaniques les plus utilisés, dresse la liste des 
paramètres de contrôle et propose une analyse dimensionnelle adaptée aux sollicitations les plus 
courantes. L'influence de ces paramètres est ensuite discutée dans différentes conditions : dé- 
formation quasi-statique de matériaux granulaires solides (§ [3]) ; écoulements granulaires denses 
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2 Simulation de matériaux modèles : paramètres de contrôle 
2.1 Géométrie 

Nous traitons presque exclusivement ici d'assemblages de particules sphériques à trois 
dimensions (3D), voire circulaires en deux dimensions (2D), avec une faible étendue granulomé- 
trique. Ces modèles font l'objet de l'immense majorité des études numériques par simulations 
discrètes. Toutefois, la forme et la distribution de tailles des grains jouent un rôle déterminant 
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dans leurs propriétés mécaniques, souvent bien davantage que certaines caractéristiques micro- 
mécaniques des contacts. C'est surtout l'expérience qui a révélé l'importance de ces facteurs 
géométriques. Quelques observations issues de travaux numériques seront néanmoins évoquées. 

D'autre part, quelle que soit la forme des grains, la géométrie des assemblages qu'ils 
forment joue un rôle central dans toutes les études micromécaniques. Ses propriétés telles que la 
densité, la connectivité, l'anisotropie éventuelle, la distribution des interstices entre grains voisins 
sont intimement couplées au comportement mécanique du matériau, qui ne résulte pratiquement 
jamais - c'est là une difficulté essentielle - d'une moyenne simple du comportement mécanique 
du contact. Dans une situation d'écoulement stationnaire (§EJ), l a microstructure détermine la 
cinématique du système autant qu'elle en est le produit. Dans un assemblage granulaire solide 
(§E]), le processus de mise en place (objet du chapitre 8) conditionne initialement la microstruc- 
ture qui influe ensuite fortement sur le comportement. Il est important, avant d'introduire les 
nombreux paramètres liés aux lois de contact, de rappeler que les propriétés mécaniques des ma- 
tériaux granulaires dépendent, autant que des propriétés des contacts entre grains, de la manière 
dont ceux-ci sont assemblés. 

2.2 Lois de contact 

Les ingrédients des lois de contact mises en œuvre en simulation numérique résultent en 
principe d'une modélisation micromécanique [T], éventuellement simplifiée et complétée par des 
hypothèses physiquement plausibles, dans laquelle certains choix sont guidés par un souci de 
faisabilité des calculs. Il convient bien sûr de s'assurer que les propriétés auxquelles on s'intéresse 
ne sont pas affectées par un choix trop arbitraire, voire qui conduirait à une violation de propriétés 
physiques fondamentales comme l'objectivité ou les lois de conservation. Le plus souvent, des 
compromis sont adoptés entre la conformité du modèle à ce qui est connu des lois de contacts 
et la simplicité des calculs, en exploitant les connaissances disponibles quant à la sensibilité 
des propriétés visées aux ingrédients micromécaniques. Dans tous les cas, il est nécessaire de 
s'appuyer sur des études paramétriques (pour lesquelles l'analyse dimensionnelle est une aide 
précieuse) . 

Nous ne donnons ici qu'une description abrégée de ces lois de contact, destinée à faire 
apparaître la liste des paramètres de contrôle influant sur les résultats des simulations. Nous 
renvoyons aux références bibliographiques, parmi lesquelles [T] est un traité de la mécanique du 
contact et 13] contiennent des introductions générales adaptées à la simulation des assemblages 
de grains, pour une présentation plus complète. 

2.2.1 Frottement 

Nous nous limitons dans un premier temps aux cas de contacts quasi-ponctuels : l'étendue 
de la région de contact est supposée faible devant le diamètre d des grains. Au point de contact 
entre les grains i et j, ceci entraîne que les efforts de contact se résument à une force F, que 
l'on décompose en une partie normale, Fn (scalaire) et une partie tangentielle F y (un vecteur 



En l'absence de cohésion, la force normale est répulsive : on a Fjy > 0. En général on 
admet l'existence d'un coefficient de frottement \i qui intervient dans la condition de Coulomb : 



2.2.2 Élasticité 

L'élasticité du contact est bien connue dans le cas d'objets réguliers sans arêtes pQ. La loi 
de Hertz relie alors, pour des grains sphériques de diamètre d constitués d'un matériau élastique 



en 3D). 
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de module d'Young E et de coefficient de Poisson u, la force normale i*V à la déflexion normale 
h des surfaces en contact (« interpénétration » apparente) : 



La variation de la partie élastique de la réaction tangentielle Ft avec le déplacement relatif 
tangentiel ô [3], du fait du critère de Coulomb écrit au niveau local du vecteur-contrainte dans 
la petite surface de contact |T], ne peut en principe s'écrire que de façon incrémentale, et donne 
des lois fort compliquées [5] . On est souvent amené à les simplifier jÏÏJ 18] , en gardant la raideur 
tangentielle Kx(h), qui apparaît dans : 

Fd 1 / 2 

dF T = h^dô = K T {h)dl (3) 

indépendante de ô tant que la condition de Coulomb ([!]) reste satisfaite. Dans le cas contraire, 
une fois incrémentée Ft dans le calcul, il faut la projeter sur le cône de Coulomb. De plus il faut 
tenir compte des variations simultanées de Fpf et de Fy, et transporter les forces de contact avec 
les mouvements d'ensemble et relatifs (y compris roulements et pivotements) des deux objets. 
Pour ce faire on ne dispose pas en fait de loi établie pour tous les cas, et on doit se contenter 
de règles praticables qui ne violent pas l'objectivité (discutée dans [9]), ni le second principe 
de la thermodynamique (il ne faut pas créer de l'énergie élastique dans un cycle, ce point est 
étudié dans [3]). Une mise en œuvre des lois de Hertz-Mindlin simplifiées qui satisfait à 

ces conditions est décrite dans [7]. 

Il arrive aussi souvent que l'on ne cherche pas à décrire précisément l'élasticité du contact, 
et que l'on choisisse plus simplement une loi linéaire unilatérale : 

F N = K N h, dF T = K T dô, (4) 

avec la condition de Coulomb ([]]). Les raideurs normale (ifjv) et tangentielle (Kt) peuvent être 
choisies pour respecter l'ordre de grandeur des déformations élastiques dans les contacts. 

Dans les relations (|3j) ou (J4j) , ô représente en fait la partie élastique du déplacement relatif 
tangentiel, qui ne varie que lorsque la force de contact est intérieure au cône de Coulomb défini 
par ([T]), ou bien pénètre à l'intérieur à partir du bord |10[ 111] . 

Parfois, on considère simplement les forces élastiques comme une manière commode d'as- 
surer l'impénétrabilité des grains et la mobilisation progressive du frottement. Une telle attitude 
implique que l'on se place implicitement dans la limite des grains rigides, sur laquelle nous re- 
viendrons. 



2.2.3 Forces visqueuses et coefficients de restitution 

La partie visqueuse de la force de contact est le plus souvent prise linéaire dans la vitesse 
relative : 

Fn = *n^, F T = a T -, (5) 

En général les coefficients d'amortissement sont assez mal connus et ils sont rarement identifiés 
à partir d'un modèle physique. (Une approche par la mécanique du contact entre solides vis- 
coélastiques a cependant été développée, que l'on trouvera présentée dans la réf. |12] et dans 
bien d'autres écrits du même groupe d'auteurs. Ce modèle a été exploité dans l'étude des « gaz 
granulaires »). 
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La modélisation d'un choc binaire frontal avec les relations (|4]) et ([5]) définit un problème 
d'oscillateur harmonique amorti, et fixe la durée r c de la collision entre deux grains de masse m : 



7T 



2K N /a N yy 1/2 

m V m J 



(6) 



Le pas de temps en dynamique moléculaire est choisi comme une petite fraction de yJm/K^. Si 
a m dépasse la valeur critique a c N = \J2mK^ (pour laquelle r c — > oo dans ©), le mouvement 
relatif lors du choc est suramorti. Comme il y a alors deux temps caractéristiques dans l'évolution 
de la déflexion normale dans le contact, dont le rapport augmente très rapidement avec £ = 
«Ar/a^r, c'est une situation numériquement défavorable, que la commodité des calculs conduit à 
proscrire. Aussi, lorsque l'on souhaite, afin d'approcher rapidement des états d'équilibre que l'on 
veut étudier, dissiper efficacement l'énergie, on tend à choisir des coefficients d'amortissement 
proches du niveau critique, mais inférieurs (£ < 1). Les lois élastique (j!]) et visqueuse ([5j) ci-dessus 
se traduisent, dans une collision isolée, par un coefficient de restitution normal e, fonction de Ç : 



exp 



TTC 



(7) 



Dans le cas hertzien, on peut se référer à la raideur linéaire tangente, K^{K), qui varie comme 
h 1 / 2 d' après (j2|), et appliquer les mêmes considérations pour le choix du pas de temps (alors 
variable selon le niveau de déflexion atteint) et de l'amortissement que dans le cas linéaire. On 
montre que le choix d'un amortissement égal, à chaque instant, à une fraction constante £ de 
la valeur critique d'un modèle linéaire avec raideur K^{h) donne un coefficient de restitution 
normal indépendant de la vitesse relative dans une collision binaire. 

Le choix du coefficient d'amortissement tangentiel ay est guidé par des considérations 
similaires à celui de (à noter que sa valeur critique fait intervenir les moments d'inertie 
des grains en contact). Sa relation avec un coefficient de restitution tangentiel est un peu plus 
complexe en raison de la fréquence différente des oscillations de h et de ô dans un choc et de la 
saturation possible du critère de Coulomb. 

Après avoir introduit des termes visqueux dans la force de contact, on doit se demander 
si on les inclut dans les composantes Fpj et Fr qui doivent satisfaire l'inégalité de Coulomb (JTJ). 
Deux attitudes sont possibles, et les deux choix ont été adoptés dans la littérature. On peut 
considérer que l'inégalité de Coulomb ne concerne que les composantes élastiques de la force 
de contact. Tout se passe alors comme si les forces visqueuses étaient dues à la présence d'un 
fluide dans la région du contact (mais sans que l'on tienne compte de la géométrie précise de son 
écoulement). La force normale totale peut alors être attractive lorsque les deux grains s'éloignent 
assez rapidement. L'autre option, qui peut sembler plus logique si on attribue les forces visqueuses 
à la dissipation interne au matériau qui constitue les grains, consiste à ajouter les forces visqueuses 
aux forces élastiques avant d'imposer l'inégalité (JTJ), auquel cas il convient d'interdire les forces 
attractives en limitant au besoin l'intensité de la force visqueuse. 



2.2.4 Autres mécanismes de dissipation 

La dissipation d'énergie dans les contacts peut aussi résulter d'un comportement plastique 
du matériau qui constitue les grains. Des modélisations de contacts élastoplastiques pour les 
objets de forme régulière (tels que les billes métalliques) sont présentées dans |T], et des lois 
très simplifiées sont mises en œuvre dans certaines simulations (ainsi dans la référence |Ï3], ou 
dans [14J où un modèle est proposé avec plasticité et cohésion). Nous ne discuterons pas des 
paramètres de ces modèles. 

D'autres mécanismes de dissipation sont utilisés comme de simples astuces nu-mériques 
destinées à favoriser l'approche d'états d'équilibre, sans grand souci de vraisemblance physique 
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(même dans la forme des lois choisies), voire en ignorant certaines lois de base. Ainsi, une force 
visqueuse est parfois mise en œuvre qui s'oppose simplement à la vitesse de chaque particule 
(comme pour un objet solide isolé baignant dans un fluide visqueux de vitesse nulle à l'infini). 
Un autre mécanisme (proposé notamment dans certains logiciels commerciaux) consiste à exercer 
sur chaque grain une force proportionnelle à son accélération, avec un coefficient de signe opposé 
à celui du produit scalaire de la force avec l'accélération. Tout se passe donc alors comme si la 
masse de chaque objet était différente dans les phases d'accélération ou de ralentissement. Ces 
procédés violent la conservation de la quantité de mouvement, et font jouer un rôle particulier à 
une vitesse de référence. On ne doit donc les utiliser qu'avec précaution, dans les cas où le bilan 
de quantité de mouvement ne joue pas de rôle important. Il faut bien sûr éviter de les mettre en 
œuvre pour modéliser un écoulement ou la mise en place d'un matériau granulaire sous gravité. 

2.2.5 Cohésion 

C'est l'énergie par unité d'aire 7 des interfaces entre les grains solides et le milieu ambiant 
qui est à l'origine de l'adhésion aux contacts, qu'elle soit ou non transmise par un ménisque 
liquide. La force attractive maximale est toujours d'ordre jd. En présence d'un pont liquide, il 
est possible d'utiliser des modèles quantitatifs pour les assemblages de billes |!5[ I16j . La portée 
Dq des forces attractives à distance est alors déterminée par la distance de rupture, de l'ordre du 
volume du ménisque élevé à la puissance 1/3. La prise en compte de l'adhésion capillaire dans la 
simulation des matériaux granulaires fait l'objet du chapitre 12. Dans le cas de grains solides secs, 
la force adhésive est sensible aux irrégularités de leur surface, et sa portée (celle des forces de van 
der Waals) est de quelques nanomètres. C'est pourquoi elle se manifeste seulement en pratique 
pour les particules de petite taille, typiquement de l'ordre du micron ou de la dizaine de microns. 
Pour des grains de sable, ou des billes dont le diamètre est de l'ordre de 100 //m, seules quelques 
aspérités peuvent être sensibles à une telle attraction, dont l'effet se trouve ainsi fortement réduit, 
alors que les forces capillaires en présence d'un liquide ne sont pas sensibles à ces irrégularités 
géométriques et entrent en jeu pour des grains plus gros. Il existe sur la modélisation des effet 
simultanés de l'adhésion et de la déformation élastique au contact une littérature fournie |17j . 
mais l'étude des assemblages granulaires cohésifs par simulation numérique (exposée au chapitre 
11) en a pour l'essentiel retenu des versions assez simplifiées, car les approches classiques restent 
incomplètes (on ne sait pas toujours bien prendre en compte les efforts tangentiels), s'appliquent 
difficilement à des grains à surface irrégulière et font appel à des informations détaillées sur 
les matériaux et les états de surface qui ne sont pas toujours disponibles en pratique. Les lois 
simples mises en œuvre dans les exemples simulations rapportés aux § [3] et H] consistent à ajouter 
un terme attractif simple |18[ 119]. à la répulsion élastique exprimée par introduisant ainsi les 
deux paramètres importants communs à tous les modèles de cohésion : une traction maximale Fq 
que peut supporter un contact, et la portée Dq de l'attraction mutuelle entre surfaces solides. A 
l'équilibre, en l'absence d'action extérieure, une paire isolée de grains est maintenue en contact 
avec une déflexion 

h ~ F /K N , (8) 

pour laquelle la force de répulsion élastique et la force d'adhésion se compensent. 

En présence d'adhésion, il est important d'appliquer l'inégalité de Coulomb ([!]) en excluant 
de Fff au second membre le terme attractif. Ainsi, dans un modèle simple où ce terme est 
constant, égal à — Fq pour les grains en contact |19j . on impose la condition ||Fj-|| < (j,(Fn +Fq). 
En effet, il faut considérer que les surfaces sont pressées l'une contre l'autre par l'attraction des 
deux grains, qui agit de façon analogue à une force extérieure. Il en résulte qu'en l'absence de force 
extérieure, le contact entre deux grains collés avec la force Fq est alors capable de transmettre 
une force tangentielle maximale égale à [iFq, alors que la force normale totale est nulle. 
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2.2.6 Résistance au roulement 

Si on prend en compte l'extension finie des régions de contact, qui ne se réduisent pas à des 
points, soit en raison de la déformablité des grains solides, soit de leur forme irrégulière et non 
convexe, on est amené à autoriser les contacts à transmettre un moment T, qui travaille dans les 
mouvements relatifs de roulements sans glissement (différence de taux de rotation autour d'un 
axe tangentiel) ou de pivotement sans glissement (différence de taux de rotation autour de l'axe 
normal). Nous évoquerons quelques exemples de l'effet d'un tel ingrédient micromécanique, encore 
relativement peu usité, dans des exemples bidimensionnels (pas de pivotement). Pour prendre 
en compte ce phénomène les modèles retenus en simulation sont restés fondés sur des approches 
simplificatrices |20| . Le paramètre essentiel est le coefficient de frottement de roulement, fj,R, 
tel que le moment maximal soit [irFn. Ainsi défini, //r a la dimension d'une longueur, et c'est 
pourquoi on rencontre aussi la définition dans laquelle c'est fj,R divisé par un rayon de grain R 
qui est appelé coefficient de frottement. La réponse d'un contact sollicité en roulement, avant que 
T atteigne sa valeur limite, peut être modélisée comme purement élastique, et on obtient alors 
un modèle analogue à celui de la loi de contact tangentielle |21[ [T9] élastoplastique habituelle. 
Dans le cadre d'un modèle de grains rigides, on peut aussi étendre l'approche de la dynamique 
des contacts au cas de la résistance au roulement [22J. 

Si le diamètre de la région du contact est l alors doit être de l'ordre du diamètre l 
de la région de contact (distance entre aspérités), tandis que la raideur en rotation relative doit 
être d'ordre K^l 2 . En présence de cohésion, si on écrit, comme pour la condition de Coulomb 
ordinaire, l'inégalité T < [J-rF^ avec la seule composante répulsive élastique de la force normale, 
on obtient des contacts capables de transmettre un moment d'ordre hrFq même si la force 
normale totale est nulle. Un faible coefficient de frottement de roulement peut donc avoir une 
influence importante dans le cas cohésif |19] , 

2.2.7 Cas des grains rigides 

La liste des paramètres que nous avons introduits jusqu'ici se trouve plus réduite dans 
le cas d'une modélisation (comme en dynamique des contacts, cf. le chapitre 3) avec des grains 
parfaitement rigides. Ainsi toutes les raideurs et les échelles de longueur associées à une déflexion 
élastique disparaissent. Restent le coefficient de frottement, éventuellement le frottement de rou- 
lement, et les coefficients de restitution, ainsi que les données ayant le sens d'une force, comme 
la traction maximale Fq, et éventuellement la portée -Do des forces attractives. 

Toutefois, comme la condition de stricte impénétrabilité n'est forcément satisfaite qu'ap- 
proximativement dans la pratique, la méthode introduit des niveaux de tolérance. Elle fait éga- 
lement appel au choix d'un pas de temps. Nous renvoyons au chapitre 3 pour le rôle de ces 
paramètres liés à la mise en œuvre numérique de la dynamique des contacts. 

2.3 Analyse dimensionnelle 

L'analyse dimensionnelle, fondée sur l'indépendance des propriétés physiques vis-à-vis du 
choix des unités fondamentales de temps, de longueur et de masse, nous assure que les propriétés 
du matériau, si elles sont exprimées par des grandeurs sans dimension (rapport de contraintes 
fonction d'une déformation par exemple), ne dépendent que des combinaisons sans dimension 
des paramètres du problème. 

Dans un essai mécanique, les paramètres, en plus de ceux du matériau, comprendront le 
plus souvent une pression P (contrainte de confinement) et un certain taux de déformation ë. 

Bien entendu, on peut construire de bien des façons différentes une liste de nombres sans 
dimension à partir des données. Le choix que nous proposons est guidé par les considérations sui- 
vantes. Ces nombres, d'abord, expriment, comme il est accoutumé, le rapport de deux grandeurs 
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physiques, par exemple à partir des échelles de temps : le temps de collision r c (qui est compa- 
rable au temps de propagation de la quantité de mouvement à travers un grain), le temps de 
cisaillement r s = 1/è et le temps inertiel T{ = y/m/P (le temps de déplacement caractéristique 
d'un grain de masse m soumis à la pression P). En présence de modèles de matériaux diffé- 
rents, on peut donner des définitions assez similaires des paramètres de contrôle adimensionnés 
qui permettent de situer dans l'espace des paramètres des régions où les comportements seront 
comparables. Il convient ensuite, afin de pouvoir faire appel le plus souvent possible aux mêmes 
définitions pour les modèles de grains rigides ou déformables, d'isoler les paramètres de raideur 
en les faisant intervenir dans le moins de combinaisons possibles. Enfin, il est préférable que les 
propriétés du système étudié possèdent des limites bien définies pour les valeurs extrêmes, très 
faibles ou très élevées, des nombres choisis. 



2.3.1 Grandeurs liées au matériau 

Le comportement des contacts fait directement intervenir des grandeurs sans dimension : 
le coefficient de frottement fj,, e (ou £) pour la dissipation visqueuse ou la restitution normale, 
ex (ou (y), son homologue pour la force tangentielle, et éventuellement le niveau de frottement 
de roulement Hr/cL. 



2.3.2 Raideurs adimensionnées 

Les raideurs normale et tangentielle se comparent entre elles, ce qui donne le paramètre 

Kt/Kn pour l'élasticité linéaire dU). Pour le modèle simplifié ©-([SI de l'élasticité de Hertz- 

2 - 2v 

Mindlin, ce rapport de raideurs est indépendant du la force et vaut — , qui varie entre 2/3 

et 1 pour les matériaux ordinaires dont le coefficient de Poisson v varie entre 1/2 et zéro. Il s'agit 
donc d'un paramètre dont les variations sont très restreintes quand on cherche à reproduire les 
propriétés des matériaux habituels. La comparaison de la raideur normale Kn avec le niveau de 
contrainte P, permet de définir un niveau de raideur k, tel que l'on ait pour la déflexion moyenne 
h dans les contacts : 

- OC K~ . (9) 

La force normale moyenne (Fn) étant proportionnelle à Pd D ~ l en dimension D égale à 2 ou 3, 
on prend 

K N 



(cas linéaire, D=2 ou 3) 
(cas hertzien, D=3) 



Pd D - 2 

E y/s (Jf) 

P(l - V 2 ) ) 

En fait, on sait relier le rapport (Fn) / PdP~ x à la densité et au nombre de coordination |16 | [T9 l [7] 
et on montre alors que dans les conditions habituelles (faible étendue granulométrique) on a un 
coefficient de proportionnalité proche de 1 dans (j5J). 

La limite des grains rigides est celle oùk4 +oo. Dans cette limite, si elle est bien définie, 
les propriétés du matériau ne doivent plus dépendre de k, et les petites déflexions aux contacts, 
d'ordre k -1 , ne modifient pas sensiblement la structure de l'assemblage granulaire par rapport à 
celle d'un assemblage de grains indéformables. 
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2.3.3 Nombre d'inertie 



Le rapport des échelles de temps imposées par la sollicitation et associées à la dynamique 
des grains de masse m soumis à des forces d'ordre Pd D ~ l définit le nombre d'inertie, 



m , . 

I = ^Pd^- W 

(Une définition très proche, dont l'écriture est identique à 2 ou 3 dimensions, consiste à prendre 
I = èd^J^, où p m est la masse volumique du matériau constituant les grains). La limite de 
l'évolution quasi-statique, dans laquelle on s'attend à ce que le système reste proche d'un état 
d'équilibre à chaque instant, est la limite où / — > 0. 

2.3.4 Paramètres adimensionnés liés à la cohésion 

L'effet principal des forces attractives est d'introduire une échelle de force, Fq, qui se 
compare aux forces liées à la contrainte de confinement. On définit ainsi un nombre de cohésion 
r\ qui caractérise l'intensité des forces adhésives relativement aux forces de confinement, ou bien 
une pression réduite F* : 

Fo * 1 d D - 1 P , , 

1-5*3?. p * = r— ' (12) 

Lorsque P* 3> 1, les forces dues au confinement dominent et les effets de la cohésion sont né- 
gligeables. La structure de l'assemblage est similaire à celle d'un système sans cohésion. Dans le 
cas contraire, P* <C 1, la cohésion peut stabiliser des structures très différentes, en particulier 
beaucoup plus lâches. De plus, c'est alors la cohésion qui détermine le niveau typique des dé- 
flexions dans les contacts, qui est d'ordre Hq définie en (|8]), et des forces normales, alors d'ordre 
Fq (positives ou négatives). Lorsque la cohésion domine, plutôt que d'évaluer l'effet de l'inertie 
avec le nombre / (|lip . on peut définir 



md , . 

'-te < 13 > 

où la force typique PdP~ x est remplacée par Fq. 

L'autre nombre sans dimension introduit par les modèles de cohésion est le rapport Do/d, 
pour lequel on ne rapportera pas ici de résultat d'étude paramétrique. 



2.3.5 Bilan et remarques. 

Il faut donc étudier, plutôt que l'influence des 8 grandeurs dimensionnées m, d, Kjy (ou 

E), Kt, è, P, Fq, Dq sur le comportement mécanique, celui des 5 paramètres de contrôle sans 

dimension k, I, P* (ou rj), Kt/Kn, Do/d, auxquels il faut adjoindre les coefficients de frottement 

et de restitution. Si l'influence des rapports Kt/Kn et Do/d se révèle assez secondaire, les trois 

nombres k, I, et P* jouent un rôle central, ainsi que les coefficients de frottement, \x et fiFi/d, 

comme il est montré dans les résultats rappelés dans les § [3] et [4] ci-dessous. Si on s'intéresse 

au comportement mécanique macroscopique, on doit provilégier l'étude des systèmes homogènes 

dans la limite où le nombre N de grains tend vers l'infini, et N peut être considéré comme un 

autre paramètre à faire varier. Comme la durée des calculs est proportionnel à N et au nombre de 

pas de temps en dynamique moléculaire, le coût de la simulation d'un intervalle de déformation 

N 

donné lorsque le taux é est imposé sous un niveau de contrainte P varie comme — y/ë. D'où 

l'importance pratique, si on s'intéresse à la limite quasi-statique, et à la limite des grains rigides, 
de savoir quelle valeur de I sera assez faible et quelle valeur de k sera assez élevée. 
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3 Déformation de matériaux granulaires solides 



3.1 Assemblage, confinement, compression 

Les matériaux granulaires à l'état solide sont le plus souvent étudiés sous un certain niveau 
de confinement. Celui-ci varie typiquement de quelques kPa à quelques MPa dans les applications 
à la mécanique des sols, ce qui correspond au poids d'une couche de sable dont l'épaisseur va 
de quelques centimètres à plusieurs centaines de mètres. Pour des billes de verre, en prenant 
E = 70 GPa et v = 0.3 dans J2J), K, d'après {TU]), décroît de 181000 pour P = 1 kPa à 1810 pour 
P = 1 MPa. 

L'état d'un matériau sous une pression de confinement P dépend en général à la fois 
du processus d'assemblage et de l'histoire de sa compression. Les méthodes numériques d'assem- 
blage, dans lesquelles un matériau solide est constitué, soit par condensation d'un gaz granulaire, 
soit par un algorithme géométrique, font l'objet des chapitres 7 et 8 du présent ouvrage. Les mé- 
thodes mécaniques, dans lesquelles on simule effectivement la trajectoire des grains pendant 
l'assemblage, font apparaître d'autres nombres sans dimension, et produisent des résultats qui 
peuvent dépendre des paramètres qui auront par la suite peu d'influence (ou aucune influence) sur 
le comportement quasi-statique. Ainsi la densité des échantillons obtenus par pluviation contrô- 
lée |23| est sensible aux coefficients de restitution normal et tangentiel, et la connectivité des 
structures lâches obtenues par un processus d'agrégation balistique isotrope de grains cohésifs 
dépend de la vitesse moyenne quadratique initiale des particules, qu'il faut comparer à la vi- 
tesse d'éloignement minimal pour échapper à l'attraction d'une particule voisine, qui est d'ordre 
\J D^F^jm |19J . L'assemblage par compression homogène à partir d'une configuration lâche et 
sans contact peut être obtenu numériquement, en particulier en recourant aux conditions aux 
limites périodiques décrites au chapitre 6. S'il est assez graduel, c'est-à-dire si la valeur maximale 
de I reste assez faible (de l'ordre de 10 -4 ), alors, avec des grains non cohésifs et un niveau de 
raideur assez élevé (k > 5000) ce procédé conduit à des configurations d'équilibre sous pression 
contrôlée qui ne dépendent que du coefficient de frottement \x (si [ir = 0) [7]. 

La compression quasi-statique isotrope de matériaux non cohésifs (c'est-à-dire la poursuite 
du processus au-delà de la formation d'un assemblage solide) entraîne une évolution modérée et 
progressive |24j de la microstructure du système. Il en est de même pour tout trajet de chargement 
tel que le rapport des contraintes reste constant, comme la compression œdométrique. Dans un 
assemblage de billes identiques, si le nombre de coordination z est élevé, c'est-à-dire proche de 
6, sous faible pression {k —> oo), il varie d'environ 5% quand k reste supérieur à 2000; s'il 
est initialement faible, proche de 4 (le plus souvent, plus près de 4,5 si on élimine les grains 
« flottants » [7J, qui représentent alors souvent plus de 10% du total), il augmente plus vite, 
approchant 4,8 pour k = 2000 (évolution facilitée par le recrutement d'une bonne moitié des 
« flottants » initiaux par la structure qui porte les forces). Pour k de l'ordre de 1000, on constate 
une augmentation considérablement plus rapide du nombre de coordination avec la pression. On 
s'attend donc à un comportement proche de la limite rigide pour k > 2000 si z est élevé, et pour 
des k un peu plus grands si z est faible. C'est la création de nouveaux contacts là où il y a avait 
sous pression plus faible de petits interstices entre grains voisins qui augmente la coordinence. On 
doit donc conclure que pour approcher de la limite des grains rigides il ne suffit pas de satisfaire 
à la condition k <C 1 que la déflexion élastique h des contacts soit faible devant le diamètre 
d des grains. Il faut en plus que ces déflexions n'entraînent que de faibles modifications de la 
distribution des interstices entre grains voisins. D'où la nécessité de satisfaire une inégalité plus 
exigeante, qui demande un niveau de raideur k minimal, dans la pratique des systèmes simples 
considérés dans ce chapitre, de l'ordre de quelques milliers. 

La compression isotrope entraîne très peu de déformation irréversible avec des grains sans 
cohésion [24j, en dépit de certaines évolutions internes (une coordinence z initialement élevée 
peut diminuer dans un cycle de compression). La situation est tout autre avec des grains co- 
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FIGURE 1 - Aspect de la configuration pour P* = 10 2 , à gauche et pour P* = 13 (et k — 7500), 
à droite, d'un échantillon 2D de disques cohésifs [25J. 

hésifs, car une augmentation de pression modifie le rapport rj = 1/P*, défini en f j 1 2 [) . entre la 
cohésion, capable de stabiliser des structures lâches, et la pression de confinement. L'application 
d'un confinement modifie alors considérablement l'état du système, la compression sous charge- 
ment isotrope (ou sur un trajet de chargement tel que le rapport des contraintes reste constant) 
conduit à une déformation plastique très importante et définit la courbe de consolidation (fi- 
gure [2]). Cette évolution de 1/$ (ou de l'indice des vides) avec le logarithme de la pression est un 
résultat classique pour les poudres et les sols cohésifs (26], et on la mesure en laboratoire dans la 
limite quasi-statique. A faible P* , l'approche de cette limite, lors de la compression, est à évaluer 
avec le paramètre I a de fj 1 3 P |25| . La courbe de consolidation est sensible à la résistance au rou- 
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FIGURE 2 - Courbe de consolidation en compression isotrope et effet de la décompression dans 
un système modèle cohésif 2D, (a est ici le diamètre maximal d'une distribution uniforme, le 
minimum étant a/2) : variations de 1/$ (ou de l'indice des vides —1 + 1/$) en fonction de 
lnP* : une faible valeur de Hr/cl suffit à en modifier la pente dans le régime linéaire. Noter la 
faiblesse de la déformation à la décharge (partie inférieure du graphe). Résultats tirés de [25J. 

lement, même très faible, \iRjd = 10~ 2 . La résistance au roulement joue un rôle micromécanique 
important dans les structures cohésives lâches |19| 125]. car elle permet la transmission d'efforts 
par des chaînes de grains, alors dotées d'un raideur en flexion. 

3.2 Matériaux non cohésifs, déformation quasi-statique 

Nous prenons l'exemple de la simulation du comportement d'un assemblage de billes iden- 
tiques non cohésives dans une compression triaxiale de révolution, essai classique en mécanique 
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des sols |26j , souvent simulé par éléments discrets avec des matériaux modèles |27[ [28] . Le taux 
de déformation axiale e a est alors imposé, et P, dans la définition de / et de K, désigne la pres- 
sion isotrope initiale, à laquelle restent égales les deux contraintes principales 02 et 03, tandis 
que 0"i = P + q augmente avec €\ — c a . Il s'agit normalement d'une sollicitation quasi-statique. 
Alors que le nombre d'inertie dans un essai triaxial de laboratoire sur un sable prend souvent 
des valeurs de l'ordre de 10 -9 , les simulations les plus soigneuses ne descendent guère en dessous 
de / = 10~ 5 . Pour que l'essai numérique soit aussi lent que l'essai de laboratoire il faudrait 
donc de l'ordre d'une année de calcul au lieu d'une heure. Les effets dynamiques peuvent donc 
se manifester bien davantage dans le calcul numérique, et il faut s'assurer qu'ils n'influent pas 
trop sur les comportements mécaniques et les mécanismes rhéophysiques observés. La figure |3] 
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FIGURE 3 - Déviateur des contraintes q, normalisé par P (axe de gauche) et déformation volu- 
mique e v (axe de droite) fonctions de la déformation axiale dans la compression triaxiale d'as- 
semblages de billes (Hertz-Mindlin, [i = 0,3, ^jlr = 0, k = 8400). On compare : (a) pour le même 
état de départ, différentes valeurs de I (Eq. HT]) et £ (cf. §2.2.3^ ; (b) deux états initiaux, A et C, 
de même densité, mais de nombre de coordination différent. 

représente la réponse au test triaxial, exprimée selon l'usage avec le déviateur des contraintes q 
et la déformation volumique e v fonctions de la déformation axiale e a , d'échantillons numériques 
de N=4000 billes, entre l'état isotrope de départ et le maximum de déviateur. 

La figure [3^ montre que, pour le cas considéré, celui d'un système dense, un nombre 
d'inertie inférieur à 10 -3 peut suffire à approcher la limite quasi-statique de la courbe contrainte- 
déformation. Les courbes macroscopiques apparaissent alors indépendantes du niveau de dissi- 
pation visqueuse (on peut même prendre C = 0) et du nombre d'inertie /. Toutefois, au niveau 
des grains, il subsiste des écarts à l'équilibre, qui dépendent de £ et tendent à se réduire si / 
décroît. Il en résulte un certain fluage lorsqu'on arrête d'imposer la valeur de e a pour maintenir 
q constant. Le système se stabilise alors finalement dans un état d'équilibre (avec les critères 
discutés au chapitre 1), après que e a et e v ont subi, typiquement une augmentation d'ordre 10~ 3 
ou 10~ 4 dans le cas considéré ici. Ces déformations différées dépendent de divers paramètres dy- 
namiques (C, /, pilotage des conditions aux limites...) mais diminuent progressivement à mesure 
que I décroît. 

La figure [3j> illustre l'importance de la géométrie de la configuration initiale pour le com- 
portement mécanique (comme annoncé au § 12 . 1 [) . Au-delà de la seule compacité, très proche pour 
A et C de la valeur maximale <&* ~ 0, 64 pour les assemblages désordonnées de billes, on voit que 
le nombre de coordination initial z détermine la forme de la courbe de déviateur, qui est beau- 
coup plus raide pour A (z ~ 6) que pour C (z ~ 4). De plus, ces deux assemblages qui diffèrent 
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par leur coordinence initiale ne présentent pas du tout la même sensibilité au niveau de raideur 
K. Le comportement du matériau C, aux échelles de déformation considérées, est indépendant 
de k (si on l'exprime avec des grandeurs sans dimension comme q/P, e v ), comme le montre la 
figure H] Les courbes relatives à différentes valeurs de k pour le matériau A diffèrent nettement 
(figure [5k.) , mais tendent à se superposer dans un intervalle de déviateur assez étendu (figure [5b) 
si on adopte l'échelle k _1 pour les déformations, en représentant q/P et K€ v comme fonctions de 
ne a . 




FIGURE 4 - Effet de k sur les courbes q{e a ) et e v {e a ) pour le système C de la figure [3] (même 
représentation). Les états initiaux modélisent ici des billes de verre sous confinement croissant, 
P = 10, 32, 100, 320 et 1000 kPa, soit k décroissant de 39000 à 1800. 




FIGURE 5 - (a) : analogue de la figure [4] pour l'assemblage A (mêmes valeurs de P). Les flèches 
classent les 5 courbes q(e a )/P et e v (e a ) selon n croissant (ou P décroissante), (b) : les mêmes 
résultats avec les déformations à l'échelle (P/Pq) 2 ^ 3 oc k -1 (Pq = 100 kPa). 

L'explication de ces différents comportements est donnée au chapitre 1, où sont distingués 
deux types de déformations, celles qui proviennent de la déformation des contacts (type I), et 
celles qui résultent de réarrangements du réseau, suite à des instabilités microscopiques (type 
II). Les systèmes bien coordonnés, comme ici l'assemblage A, possèdent un réseau de contact 
initial capable de résister à l'application d'un important déviateur, et donc présentent une ré- 
ponse de type I, qui, bien que non élastique |29[ [30] en raison de la mobilisation du frottement 
intergranulaire, est telle que la déformation soit inversement proportionnelle au niveau de rai- 
deur. La rupture de la structure des contacts survient à l'approche du maximum de déviateur, 



12 



et on observe alors des déformations de type II. Lorsque, comme dans le cas de l'assemblage C, 
le réseau des contacts est mal coordonné et cède facilement, la déformation de type II domine 
pour des déviateurs beaucoup plus faibles, et c'est la géométrie du paquet de grains en contact 
qui régit l'amplitude des déformations, devenues insensibles à la raideur k. L'amplitude des dé- 
formations est alors liée aux épaisseurs des interstices entre grains voisins, dont la fermeture crée 
de nouveaux contacts et stabilise une nouvelle configuration du réseau après que l'incrément de 
contrainte (ou de déformation) imposé a conduit à l'instabilité et à la rupture de la précédente. 
En posant que les déformations d'un réseau stable (type I), sont inversement proportionnelles au 
niveau de raideur, nous admettons que sa géométrie n'y est pas sensible, ce qui est la définition 
de la limite des grains rigides. On voit donc que la limite des grains rigides correspond à la 
séparation d'échelle entre les déformations de type I (à l'échelle k" 1 ) et les déformations de type 
II (dont l'échelle est fixée par les interstices). 

Le domaine élastique, accessible à l'expérience, est strictement inclus au régime de défor- 
mation de type I, plus étendu. Les déformations de type I, bien que d'ordre k -1 , restent donc 
mesurables |30j . et on doit prendre la bonne valeur de k pour les modéliser correctement en 
simulation. C'est dans ce régime que la déformation peut aussi dépendre de Kt/Kn (ou du 
coefficient de Poisson v du matériau constituant les grains). Cependant, les résultats de |31| sur 
les assemblages de disques (avec un intervalle de déviateur important en régime I) montrent 
que le comportement mécanique est à peu près insensible à ce paramètre s'il dépasse 1/2. Une 
curiosité : le choix (plutôt injustifié pour des grains solides ordinaires) de Kt > Kn peut don- 
ner des matériaux à coefficient de Poisson négatif |32j . Quant aux raffinements des lois ((2J)- ([3j) 
proposés dans [5], on a pu vérifier [30j qu'ils n'influent que très peu sur les propriétés élastiques 
des assemblages de sphères. 

Quand on prépare un échantillon numérique dense en supprimant le frottement dans la 
phase d'assemblage (chapitre 8), on obtient une configuration de coordinence élevée, analogue 
à l'état A pris en exemple ici (figure |5]). Dans les résultats expérimentaux sur les sables, la 
croissance de q avec c a est très souvent beaucoup plus lente que celle que l'on obtient alors dans 
la simulation. Il est alors tentant d'attribuer aux raideurs choisies pour les calculs une valeur 
trop faible (par exemple K = 100). Ce faisant, on décale effectivement le maximum de déviateur 
vers des déformations plus grandes (figure |5^), mais on commet l'erreur d'adopter un modèle 
physiquement non pertinent pour les déformations. A cet égard il est préférable de choisir l'état 
initial noté C dans les exemples précédents (figure |3]b) . 

Les résultats rappelés ici (figure |3b) indiquent aussi que le rapport q ma . x /P du maximum de 
déviateur (le « pic ») à la contrainte latérale ne dépend pas de K (ici > 1800) et est essentiellement 
sensible, pour un matériau donné, à la compacité initiale $ de l'échantillon (avec ici <&c — = 
<&*), conformément aux résultats de la mécanique des sols [26]. L'angle de frottement interne <p 
est fonction croissante du coefficient de frottement intergranulaire, comme l'indique la table HJ 
qui fournit également les valeurs de la dilatance T> lorsque q = ç/ max - Les données de la table [T] 
sont donc définies par : 



Ces résultats sont relatifs à la compacité maximale <3?* des configurations désordonnées de sphères 
d'une seule taille. Il est traditionnel de relier (p à cp c et à ï>, le frottement intergranulaire et la 
dilatance apparaissant comme les origines du frottement interne |33] (voir aussi [34]). Cependant, 
il faut noter (table [T]) que les assemblages de grains non frottants sont dotés de frottement 
interne mais dépourvus de dilatance [35] . C'est la microstructure des assemblages granulaires, 
qui, en développant, sous l'effet de la déformation, une certaine anisotropie [36], permet au 
matériau de supporter des contraintes déviatoires croissantes |37 [ I38 [ l39 | 140] . Il s'agit en général 
d'un phénomène de réarrangement géométrique, associé à des déformations de type IL On a pu 
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26.6 


45 




5.5 ±0.5 


24 ± 1 


26 ± 1 


28 ± 1 




0±0.01 


0.52 ±0.01 


0.60 ±0.01 


0.71 ±0.01 



TABLE 1 - Angle de frottement interne 99 et dilatance T> pour q = g max , fonctions de [i = tgip c , 
pour des assemblages de sphères monodisperses (fin = 0). Valeurs issues de simulations de 
compressions triaxiales à partir d'états initiaux isotropes, avec $ ~ <£>* |23| . 



noter que l'introduction d'un faible frottement de roulement pouvait sensiblement augmenter le 
niveau de déviateur maximal |41[ 134] . Avec des particules ellipsoïdales plutôt que sphériques, le 
frottement interne augmente aussi très notablement [42J. 



3.3 L'état critique 



Les simulations confirment |27[ 138) 1^3] la théorie de l'état critique |26] : pour de grandes 
déformations, l'état interne du matériau s'approche d'un certain attracteur, acquiert une certaine 
structure d'écoulement, indépendamment de l'état initial. Une fois cet état critique atteint, la dé- 
formation plastique peut augmenter indéfiniment, de façon monotone, à déviateur de contraintes 
et densité constants. L'état critique correspond à la limite quasi-statique des écoulements granu- 
laires, objets du § [4] Lorsque la déformation se localise dans une ou plusieurs bandes de cisaille- 
ment, le matériau se trouve dans l'état critique à l'intérieur de ces bandes |44| . Relativement bien 
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FIGURE 6 - Frottement interne à l'état critique, /i* (a), et rapport $ c /$* de la compacité critique 
à la compacité maximale des assemblages aléatoires (b) (<£»* ~ 0,639 en 3D, et $* ~ 0,85 en 
2D). Les résultats proviennent de différentes études numérique publiées, certaines pour des grains 
d'une seule taille, d'autres introduisant une certaine polydispersité (diamètres uniformément 
distribués entre (1 — x)d et (1 + x)d , pour des systèmes 2D (symboles pleins) : |22) (carrés, 
x = 0, 11) ; [34] (triangles, x = 0, 6) ; et en 3D (symboles ouverts) : [3] et [35] (cercles, x = 0) ; 
[27] (triangle, x ~ 0, 6) ; 03] (carrés, x = 0) ; 05] (étoiles, x = 0). 



documenté en simulation en raison de son indépendance vis-à-vis des questions délicates d'as- 
semblage et de conditions initiales, l'état critique est indifféremment étudié par des méthodes 
numériques qui prennent ou non en compte la déformabilité des grains, et est insensible au niveau 
de raideur k s'il est assez élevé (supérieur à quelques milliers). Il est caractérisé par le frottement 
interne et la densité critiques, qui ne dépendent alors que de facteurs géométriques et du coeffi- 
cient de frottement \x (si /zr = 0). La figure [6] donne le frottement interne et la densité à l'état 
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critique0 en fonction de /i pour des systèmes de disques (2D) ou de sphères sans résistance au 
roulement. 

L'état critique est très sensible au frottement de roulement |22] , La référence |46j en donne 
une étude pour un matériau (2D) de grande polydispersité. Il y est montré en particulier que 
l'angle de frottement interne ne dépend pas sensiblement de la granulométrie. 



4 Écoulements 

On vient de voir que l'état critique des assemblages granulaires non cohésifs, dans lequel le 
matériau se déforme continuellement, dans la limite quasi-statique, en maintenant une structure 
interne constante, dépend exclusivement, parmi les différents paramètres de contrôle introduits 
au § 13.31 du coefficient de frottement (pour des k assez grands). Nous nous proposons maintenant 
d'étudier les écoulements stationnaires uniformes |47], dans lesquels l'état du système a ceci de 
commun avec l'état critique qu'il est maintenu constant, indépendant de l'état initial, et ceci de 
différent que l'on n'est plus dans la limite des faibles nombres d'inertie. 

La loi de comportement va alors s'exprimer à travers les dépendances des variables ca- 
ractérisant l'état stationnaire du système en fonction des paramètres de contrôle sans dimension 
précédemment introduits. On commence par le cas de matériaux granulaires non cohésifs dans la 
géométrie la plus simple, celle du cisaillement plan sans gravité où la distribution des contraintes 
est homogène |48 [ I49 [ ILS ] 150] . On discute alors les influences respectives de l'état de cisaillement 
(à savoir le nombre inertiel I) et des caractéristiques mécaniques des contacts (à l'exclusion tou- 
tefois du frottement de roulement) sur la loi de comportement. On aborde ensuite l'influence du 
gradient de contraintes |51j . et enfin l'influence de la cohésion |!8[ 150] . à travers le nombre rj. 



4.1 Cisaillement homogène 

Dans cette configuration générique, le système occupe une cellule de longueur L et de 
hauteur H, avec des conditions aux limites périodiques dans la direction du cisaillement (voir 
chapitre 6). On impose la pression latérale Pp, plutôt que la compacitéd, de sorte que la hauteur 
H n'est pas absolument fixée. On cherche par ailleurs à imposer le taux de cisaillement é. Deux 
types de condition aux limites ont été testés dans la direction perpendiculaire au cisaillement 
(voir chapitre 6). La première méthode consiste à placer le matériau entre deux parois rugueuses 
et parallèles. Pour imposer le cisaillement, une des parois est fixe et l'autre se déplace à la 
vitesse V. La seconde méthode, dite de Lees- Edwards (on parle aussi de conditions aux limites 
bipériodiques) [54] . consiste à répéter périodiquement la cellule dans la direction y orthogonale à 
l'écoulement. Pour imposer le cisaillement, les cellules inférieure et supérieure sont animées d'une 
vitesse ±l/@. Le contrôle de la pression P est assuré en permettant la dilatation de la cellule 
selon y au cours du temps. 

On évoque dans la suite les résultats obtenus pour des systèmes bidimensionnels (les 
contraintes sont alors homogènes à une force par unité de longueur), simulés par la méthode 
de dynamique moléculaire |48[ |4"9| [TH) 150] . On en tire des enseignements généraux, vérifiés dans 

1. Les données de la figure [6] correspondent à l'état critique en cisaillement simple pour les systèmes 2D et 
les systèmes 3D de [45] . [35] et |44| . Celles de [27] et de [3] ont été obtenues en compression triaxiale. L'angle de 
frottement interne peut être différent dans ces deux états - c'est la cas en particulier pour la valeur en = |40| . 

2. On constate, comme dans d'autres géométries, que les contraintes normales sont approximativement 

égales. 

3. Ceci est approprié pour discuter de situations réelles puisque les conditions expérimentales imposent 
habituellement le niveau de contrainte plutôt que la compacité. On a montré que les simulations à pression 
contrôlée sont en très bon accord avec celles à compacité contrôlée [521 153] . 

4. La vitesse V de la paroi ou des cellules images s'adapte aux fluctuations de hauteur pour maintenir un 
taux de cisaillement constant. 
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d'autres simulations (assemblages de sphères |35j . méthode de dynamique des contacts [55J). 
Sauf mention contraire, le matériau granulaire est une collection dense de N disques dissipatifs 
de diamètre moyen d et masse moyenne m. Une petite polydispersité de ±20% est introduite 
pour empêcher la cristallisation. 

Comme on souhaite étudier la rhéologie intrinsèque du matériau et non les effets de paroi, 
les parois sont systématiquement assez rugueuses afin d'éviter une localisation du cisaillement à 
leur voisinage. Elles sont constituées de grains jointifs de mêmes caractéristiques que les grains 
cisaillés, assemblés pour former un seul solide. 

Les propriétés mécaniques des contacts sont décrites par les trois paramètres mécaniques 
indépendants introduits au § 12.21 : le coefficient de frottement /i, le coefficient de restitution dans 
les collisions binaires e et la raideur Kjy {Kt est du même ordre de grandeur que Kn [I], et 
comme il a une très petite influence sur les résultats |56[ 157] . il a été fixé à Kn/2 dans tous les 
calculs). Sauf mention contraire, les résultats rapportés ici correspondent à /i = 0, 4, e = 0, 1 et 
k = 10 4 . Ceci constitue le matériau de référence pour la suite de la discussion. La longueur L 
est choisie supérieure ou égale à 40 grains pour ne pas avoir d'effet de taille |58| . L'épaisseur H 
varie entre 20 et 100 grains, de sorte que N varie entre 1000 et 5000. 

Le système évolue vers un état de cisaillement stationnaire et uniforme, indépendant du 
mode de préparation (état lâche par dépôt aléatoire, ou dense par compaction cyclique sans 
frottement), caractérisé par des profils stationnaires de compacité, vitesse et contrainte. On 
considère que les distributions statistiques des quantités d'intérêt (structure, vitesses, forces. . .) 
sont indépendantes de la direction du cisaillement et du temps, de sorte que l'on moyenne à la 
fois dans l'espace (selon la direction du cisaillement) et dans le temps sur une durée At > IO/7. 

4.2 Loi de comportement 

4.2.1 Influence du nombre inertiel 




FIGURE 7 - Évolution du réseau de contact en cisaillement plan homogène (taux de cisaillement 
et pression imposés) en fonction du nombre inertiel / : les traits représentent les forces normales 
entre particules (d'après [59J). 

L'analyse dimensionnelle a montré que l'état de cisaillement est caractérisé par le nombre 
inertiel /, grandeur fondamentale pour décrire le comportement rhéologique des matériaux gra- 
nulaires cisaillés. Une petite valeur de / correspond à l'état critique décrit au § 13.31 Inversement, 
une grande valeur de / correspond au régime collisionel, où les grains interagissent par des colli- 
sions binaires que l'on peut considérer comme instantanées et non corrélées |60j . Entre ces deux 
régimes existe un régime de cisaillement dense pour lequel l'inertie des grains n'est pas négli- 
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geable. Le matériau est au-dessus de son seuil d'écoulement, dans un état liquide plutôt que 
gazeux. Les mouvements de grains sont fortement corrélés, et l'hypothèse de chaos moléculaire 
propre aux gaz dilués est alors en défaut. Il existe un réseau de contact percolant, extrêmement 
fluctuant à la fois dans l'espace et dans le temps. Plus le nombre inertiel augmente, moins les 
chaînes de force sont longues et nombreuses (figure [TJ. Dans les états de cisaillement stationnaire 
et uniforme, à pression et taux de cisaillement imposé, la compacité et la contrainte de cisaille- 
ment s'ajustent en réponse au nombre inertiel imposé, à travers une loi de dilatance et une loi 
de frottement qui décrivent les variations approximativement linéaires des deux grandeurs sans 
dimension, la compacité 3> et le frottement effectif /i* = S/P, en fonction du nombre inertiel / 
(Fig. E} : 

~ $max " al, 

Les paramètres & mauX , M* in , a e ^ dépendent des caractéristiques du matériau granulaire. Pour 
le matériau de référence, $ max — 0,82, /i ^ i T1 ~ 0,26, a ~ 0,37 et b ~ 1,1. Dans le régime 




I 

FIGURE 8 - (a) Loi de dilatance, (b) loi de frottement (fj, = 0,4 - e variable) - d'après 1491. 



quasi-statique (J < 10~ 2 ), le frottement effectif minimal /x^, et la compacité maximale $ max 
s'identifient au frottement interne tgcfi et à la compacité $ c dans l'état critique |38| 135], Lorsque 
/ augmente, le milieu se dilate légèrement, la proportion de collisions augmente par rapport aux 
contacts maintenus jusqu'à un régime purement collisionnel (/ > 0, 2) dans lequel la dilatance 
devient plus importante et le frottement effectif semble saturer. Ainsi, à taux de cisaillement 
imposé, l'utilisation du nombre / permet de traiter les comportements liquides et solides et le 
passage entre les deux, sans singularité. 

Avec ces deux lois on peut réécrire la loi de comportement du matériau dans le régime in- 
termédiaire sous forme viscoplastique, avec un seuil d'écoulement de Coulomb et des contraintes 
visqueuses, dont l'origine microscopique se situe dans l'organisation du réseau de contacts. Une 
formulation tensorielle de ces lois été proposée et appliquée au cas d'écoulements vraiment tri- 
dimensionnels, pour lesquels la comparaison avec les résultats expérimentaux s'est avérée très 
satisfaisante [61J. 

4.2.2 Autre choix des paramètres de contrôle 

Certains auteurs |57| ont choisi de contrôler la compacité $ plutôt que la pression, et en 
conséquence, plutôt que I et k, ils ont utilisé le couple de nombres sans dimension $ et a = 
7 / ' yj 'Kn /m = I/y/n, comme variables caractérisant l'état du matériau granulaire en cisaillement 
stationnaire uniforme, a peut être vu comme le rapport des temps de collision et de cisaillement, 
ou des vitesses de cisaillement et du son (nombre de Mach |57] ) . Les deux choix sont parfaitement 
légitimes, puisque l'analyse dimensionnelle prédit soit fi* = fi(I,n) et $ = f2(I,K), soit /i* = 
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/3($, a) et k = fn{^, a), les deux résultats étant également valables. Le choix de / et de k peut 
cependant sembler plus pratique pour plusieurs raisons. D'abord, les variations des résultats avec 
$, considéré comme paramètre de contrôle, sont extrêmement rapides. Chaque matériau possède 
une compacité critique <]? c , au-dessus de laquelle il ne s'écoule plus, sauf si les contraintes sont 
si grandes que la compression élastique des contacts compense la différence <£> — <3? c . Par contre, 
en-dessous de <Ê> C , un milieu granulaire continuellement cisaillé est libre de s'écouler avec une 
contrainte de cisaillement négligeable, sauf si la vitesse est assez grande pour conduire à une 
pression significative. On doit donc piloter <3? avec une grande précision pour observer des niveaux 
de contrainte usuels. Ceci rend les comparaisons entre différent systèmes granulaires difficiles, 
puisqu'il faudrait connaître à l'avance la valeur de la compacité critique pour chacun d'entre eux. 
De plus, la limite des grains rigides devient singulière, puisque toutes les valeurs de $ au-dessus 
de <3> c sont strictement interdites pour 7 / 0, alors que toutes les propriétés des cisaillements 
avec <ï> < <1> C sont proportionelles à une puissance de 7 dans cette limite |55j . Ainsi, aucun 
changement de régime n'est attendu par un changement de 7, ce qui semble contredire l'intuition 
sauf si l'on rappelle que la limite des grains rigides, lorsque l'on augmente le taux de cisaillement, 
finit par requérir des raideurs de contact déraisonnablement élevées, du fait de la très grande 
valeur de la pression. A l'inverse, si l'on utilise I et k comme paramètres de contrôle, aucune 
singularité n'intervient dans les résultats dans les limites / — > ou k — > 00, et des matériaux 
différents devraient manifester des comportements similaires (à défaut d'être quantitativement 
identiques) pour les mêmes valeurs de ces paramètres (ce qui définit ainsi approximativement 
des états correspondants). 

4.2.3 Influence de l'élasticité des contacts 

En le faisant varier entre 40 et 2, 5 • 10 5 , on observe que k n'a pas d'influence sur la 
loi de comportement dès qu'il dépasse 10 (ou même 10 3 ), ce qui est conforme à l'observation 
dans l'état critique. Pour les petites valeurs (40), on observe une localisation près de la paroi en 
mouvement dans les grands systèmes [48j, que l'on interprète comme un effet de la décroissance 
de la longueur de corrélation du champ de déformation quand les grains deviennent mous. On 
remarque que k = 10 4 correspond à Tj/t c = 10 2 . Puisque / = Tj/r s est plus petit qu'environ 0, 1 
dans le régime dense, t s est au moins 1000 fois plus grand que r c dans la limite des grains rigides. 
En conséquence, on peut bien parler de limite des grains rigides dans le régime dense. De fait, 
les résultats obtenus par dynamique moléculaire sont en parfait accord avec ceux obtenus par la 
méthode de dynamique des contacts, correspondant à des contacts parfaitement rigides |53j . 

La prise en compte de l'élasticité des grains conduit à tracer un diagramme des régimes 
en fonction de I/\/~k et [57], mettant en évidence un régime quasi-statique avec des effets 
de l'élasticité, un régime purement inertiel sans effet de l'élasticité, et un régime intermédiaire 
élastique-inertiel, qui ne s'avère accessible que pour des grains très déformables (k < 100). 

4.2.4 Influence du coefficient de frottement 

On a fait varier le coefficient de frottement \i entre et 0, 8. On retrouve la même forme 
qualitative des lois de dilatance et de frottement, avec des valeurs extrêmes $ max et ^ m j n dépen- 
dant sensiblement de ju, conformément aux variations des deux paramètres tg0 et <3? c de la loi de 
comportement dans l'état critique, montrées au § 13.31 La compacité reste une fonction linéaire 
décroissante de /, avec a ~ 0, 35, et le frottement effectif reste une fonction linéaire croissante 
de I, avec b ~ 1 [48J. L'approximation linéaire est cependant en défaut pour I < 0, 01 dans le 
cas des grains non frottants (voir figure [9}. La limite quasi-statique, dans de tels matériaux, est 
approchée pour des valeurs beaucoup plus faibles de / que dans le cas de grains frottants, et 
A^min es t lui-même beaucoup plus petit [35] , 
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FIGURE 9 - (a) Loi de dilatance, (b) loi de frottement (fj, = - e = 0, 1 (o) - e = 0, 9 (carrés 
noirs) - d'après 1491. 
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4.2.5 Influence du coefficient de restitution 



On a fait varier e entre 0, 1 et 0,9 |48] , Dans le cas de grains frottants, le coefficient de 
restitution e n'a pas d'influence sur la loi de comportement (figure [8J. Dans le cas de grains 
non frottants, le coefficient de restitution e n'a pas d'influence sur la loi de comportement pour 
I < 0, 1 (figure [9J. Au-delà, le comportement dépend de e : lorsque e — > 1, la chute de la 
compacité est moins marquée et le frottement effectif sature (voir aussi |57j). 

4.2.6 Limite collisionnelle 

La situation / >0,l,/i = 0,e = 0,9etr/ = correspond à la limite dense de la théorie ci- 
nétique, c'est à dire des collisions binaires quasi-élastiques (la durée moyenne des contacts devient 
égale au temps de collision |52]). Dans cette limite dense, on déduit de la théorie cinétique [62] 
la valeur du frottement effectif, indépendant de / mais dépendant de e : 



Pour e = 0,9, on trouve ^* = 0, 29 ce qui est en assez bon accord avec la valeur mesurée. 
Cette évolution vers le régime collisionnel explique ainsi la saturation du frottement effectif pour 
J > 0, 1 lorsque e = 0, 9 (figure [9b). Ainsi, comme dans l'état critique, la loi de comportement 
des écoulements denses ne dépend que de / et de ^ dans la limite des contacts rigides {k assez 
grand). Il faut atteindre le régime collisionnel pour que e joue aussi un rôle. 

4.3 Influence du gradient de contrainte 

Il est indispensable de tester si la loi de comportement identifiée en cisaillement homo- 
gène est valable dans d'autres géométries de cisaillement où la distribution des contraintes est 
hétérogène. On a ainsi considéré le cisaillement stationnaire et uniforme dans la direction du 
cisaillement, dans des géométries simples où l'on connaît la distribution des contraintes. 

Sur plan incliné, une épaisseur H de matériau s'écoule sur un socle rugueux, incliné de 
9 |58[ HH]. Dans ce cas, le coefficient de frottement effectif est constant dans la couche en écou- 
lement, égal à tg#, de sorte que la compacité et la nombre inertiel sont constants dans toute la 
couche (sauf près du socle rugueux et de la surface libre). Au coeur du matériau, on retrouve 
la loi de comportement mesurée en cisaillement homogène. Le taux de cisaillement dépend du 
coefficient de frottement [63, 56j. Il est peu sensible au coefficient de restitution, si ce dernier 
reste inférieur à 0,7 |63[ 164"! 156] , Cette indépendance s'observe aussi vis à vis du coefficient de 
restitution normal au socle sur toute la gamme des coefficients de restitution |65[ 63] ■ 

En cisaillement annulaire |51| . le matériau est confiné entre deux cylindres rugueux, inté- 
rieur de rayon Ri et extérieur de rayon R e assez grand (> 2Ri) pour que la paroi extérieure n'ait 
pas d'influence. Le cylindre intérieur tourne à une vitesse de rotation Çl. Le cylindre extérieur 
exerce une pression radiale P. L'état de cisaillement est décrit par Ri et, à travers les valeurs 
imposées de P et de fi, par la vitesse tangentielle adimensionnée de la paroi Vg = - J^y^Fj *ï u * 
généralise, à l'échelle du système global, la notion de nombre inertiel introduite en cisaillement 
homogène. Une petite valeur de Vg correspond au régime quasi-statique, tandis qu'une valeur 
élevée correspond au régime collisionnel. On utilise des conditions aux limites périodiques en 
se restreignant à un secteur angulaire O tel que L = i?j© soit supérieur à une quarantaine de 
grains. Par ailleurs, les caractéristiques du système simulé sont identiques à celles du système 
en cisaillement homogène. Toutefois, si la pression radiale s'avère approximativement constante, 
la contrainte de cisaillement décroît comme l'inverse de la distance radiale r au carré. Ceci, 
combiné à des effets d'interfaces spécifiques |66) . conduit à une localisation du cisaillement près 
du cylindre intérieur. La mesure des relations entre les grandeurs locales $(r), fi*(r) et I(r) 




(16) 
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montr que le comportement rhéologique ne dépend pas de la distribution des contraintes dans 
le régime inertiel. En revanche, lorsque (jl* devient inférieur à /J min , les lois de dilatance et de 
frottement (|15p ne sont plus valables. On observe un régime de nuage dont la compréhension 
reste un sujet d'étude. Alors qu'en cisaillement homogène, /-i min est la valeur maximale de 
supportée par le matériau granulaire avant qu'il ne commence à s'écouler quasi-statiquement, il 
est capable de s'écouler sous ce seuil, lorsque la distribution de contrainte est hétérogène. 

4.4 Influence de la cohésion 

0.01 | ► . 



a 



ni 

FIGURE 10 - Influence de I et rj sur la microstructure (d'après }50l). 

Les résultats, issus de |18| 150] . que nous rappelons ici font suite au travail entamé par 
plusieurs équipes (67] ES ESI EÔJ EU EU ES El], et visent à préciser quantitativement les lois 
rhéologiques en présence de cohésion. En dehors de l'introduction du nombre de cohésion 77, les 
caractéristiques du système simulé sont les mêmes que celles décrites au § 14.11 On considère un 
modèle simple de force normale de cohésion |75[ 169] qui prend en compte la force attractive 
maximale Fq : F a (h) = —\/ÂKjyF~oh (et Dq = 0). En fait, la forme précise de la force attractive 
F a (h) semble ne pas avoir d'influence sur le comportement rhéologique, seul compte .Fn |76[ [18] . 
La composante attractive de la force normale tend à rapprocher les grains et la déformation 
moyenne h* dépend de l'intensité de la cohésion h*{rf) = j:T~L(rj) avec 

d'ordre 1 aux faibles 77 et d'ordre 4rj aux forts 77. Le paramètre k a été fixé à 10 5 , de sorte que 
l'on reste dans la limite des contacts rigides au moins pour les faibles cohésion : h*(rj) < 10~ 4 
pour 77 < 2, 5. L'influence de la déformation des grains, pour les valeurs plus élevées de 77, est 
discutée dans [73J. Lorsque l'intensité de la cohésion n'est pas trop forte (77 < 100), le matériau 
est cisaillé de manière homogène (figure [TU]) . Les lois de dilatance et de frottement (I15p conservent 
la même forme mais les paramètres ^maxi Mmin) o et 6 dépendent de 77. La cohésion provoque 
une diminution de la compacité, mais une augmentation du nombre de coordination, et une 
homogénéisation des directions de contact. C'est le signe d'une organisation des grains en amas 
compacts séparés par du vide. La figure [TTr trace les deux paramètres $ max et a en fonction 
de l'intensité de la cohésion 77. Les deux fonctions ont une forme similaire : une importante 
diminution dès les premiers niveaux de cohésion (77 < 2) puis une diminution plus faible. La 
diminution de 0(77) jusqu'à correspond à des systèmes dont la compacité ne dépend plus de 
l'état de cisaillement /. Les deux fonctions /J m j n (î?) et b{rf) ont la même forme (figure [TTb) . Sous 
un seuil de cohésion (77 > 10), la cohésion n'affecte pas /J min ou b. Au-dessus du seuil, /-t mm (ï7) 
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et 6(77) augmentent fortement. On peut ainsi distinguer des régime de faible (rj < 10) et forte 
(77 > 10) cohésion. 




a l-i 




FIGURE 11 - (a) ^max(v) (carrés noirs) et a{rf) (o) ; (b) /^mm( r ?) (carrés noirs) et b(rj) (o) 
(d'après J59j). 



La force de frottement entre les grains est décrite par un critère de Coulomb qui considère 
uniquement la composante élastique de la force normale. Ecrit en fonction de la force normale 
totale, le critère s'écrit llF^H/i^/v < /^(l-^b/^ivl)- Le coefficient de frottement apparent en 
présence de cohésion est donc : fïH(\Fo / F^\) . Pour une faible cohésion, la force normale peut 
être beaucoup plus grande que la résistance à la traction Fn 3> Fq. Dans ce cas, la fonction 
H tend vers 1 et le coefficient de frottement apparent reste \i. Pour une forte cohésion, la force 
normale est beaucoup plus petite que la résistance à la traction : Fjy <C Fq et la fonction % 
tend vers A\Fq/F^\. Par conséquent, la cohésion augmente fortement le coefficient de frottement 
apparent entre les grains, et comme ceci a été montré pour les grains sans cohésion, ceci augmente 
la dilatance du matériau. C'est ce que l'on mesure sur la figure [TÏÏl qui compare la compacité de 
grains frottants et non frottants (/j = 0). 
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FIGURE 12 - &(r],I = 0,2) pour des grains frottants fi = 0, 4 (carrés noirs) et non-frottants 
fj, = (o) (d'après JÏÏÔf). 



5 Conclusion 

Le choix de paramètres appropriés pour étudier des comportements mécaniques de ma- 
tériaux granulaires par simulation numérique discrète dépend bien sûr du problème particulier 
auquel on s'intéresse. Il peut utilement s'appuyer sur l'analyse dimensionnelle et sur les connais- 
sances que l'on a acquises, pour des matériaux modèles simples dans différents régimes de com- 
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État critique, I < 10 -3 
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Ecoulements denses, 10 — 3 < I < 0.1 
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Régime collisionnel, J > 0.1 





N 















TABLE 2 - Influence des paramètres sans dimension sur les comportements mécaniques dans 
différents régimes de sollicitation. (0=Oui, N=Non) 

portement, de l'influence des paramètres de contrôle sans dimension choisis au § 12.31 : les plus 
importants sont le nombre d'inertie / (défini en HT]) , le niveau de raideur k (défini en HU|) . et 
le cas échéant, le nombre de cohésion r\ (défini en llip. Le choix de paramètres sans dimension 
retenu permet aussi commodément de traiter des matériaux solides comme en écoulement. La 
loi constitutive en situation d'écoulement stationnaire prend la forme d'une « équation d'état » 
qui donne l'ensemble des caractéristiques mécaniques et de microstructure en fonction d'un seul 
paramètre (/) ou de deux (/ et r\) pour les systèmes cohésifs. La table |5] résume les résultats 
rapportés dans les § [3] et |U Aux trois nombres k, I et r) (ou P*) il faut ajouter le coefficient de 
frottement intergranulaire [i (et éventuellement fin/d, le frottement de roulement), dont la valeur 
influe sur l'ensemble des comportements évoqués. On doit aussi garder à l'esprit le rôle crucial de 
la géométrie, celle des grains eux-mêmes, et celle de l'état initial pour tous les comportements de 
type solide, sauf les états critiques. Le choix des paramètres I et k permet de traiter de la limite 
quasi-statique et de celle des grains rigides, importantes en pratique. La valeur / = 10 -3 indiquée 
dans la table pour l'approche de la limite quasi-statique (comportement indépendant de I) est 
assez indicative, elle dépend bien sûr des systèmes étudiés et des configurations auxquelles on 
s'intéresse. En général, la limite quasi-statique est d'autant plus facile à approcher que le système 
est dense et bien coordonné, et demande davantage de précautions en présence de déformations 
de type II avec des réarrangements importants. De même, la limite des grains rigides peut se 
situer un peu au-delà du seuil k > 1000 indiqué. L'effet du paramètre v ou de Kt/Kn est limité 
au régime de déformations de type I, et faible si Kt/Kn > 1. Dans cette limite, il peut être 
commode de modéliser les grains comme des solides indéformables, à la manière de la méthode 
de dynamique des contacts (Chapitre 3). Cette approche peut donc s'appliquer à l'ensemble des 
cas recensés dans la table, sauf le régime de déformation de type I. Celui-ci apparaît comme un 
domaine de comportement rigide indéformable s'il est traité en dynamique des contacts. Même 
si ce point mériterait des investigations complémentaires, on s'attend à ce que l'instabilité des 
réseaux de contact qui borne l'intervalle de comportement de type I coïncide en dynamique des 
contacts et dans une modélisation prenant en compte un niveau de raideur finie. L'approche par 
des grains indéformables est alors légitime tant que l'on peut ignorer les déformations de l'ordre 
de quelques ft _1 qui surviennent dans le matériau réel. 

Les processus d'assemblage figurent dans la table, mais il conviendrait bien sûr de distin- 
guer les différentes méthodes de préparation d'échantillons solides (chapitre 8). Leur description 
précise demande d'introduire d'autres paramètres de contrôle. Comme il s'agit souvent de pro- 
cédés impliquant le passage d'un état agité et lâche à un assemblage solide en équilibre, les 
paramètres dynamiques I et e sont susceptibles d'influencer le résultat (comme dans le cas de 
la pluviation contrôlée |23j). La compression de matériaux solides dans un trajet de chargement 
conservant le rapport des contraintes, comme évoqué au § 13. 1[ est à classifier avec le comporte- 
ment solide en déformation de type I en l'absence de cohésion, et comme déformation de type 
II dans le cas cohésif (on évaluera alors, si rj > 1, les effets dynamiques avec I a ). En présence 
de cohésion, tous les comportements sont sensibles à la valeur de r], sauf bien sûr dans la limite 
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rj — > (pour la consolidation 7] < 5.10 -2 semble suffisant pour négliger la cohésion). 

Au-delà de son utilité pratique pour le choix de paramètres appropriés dans une application 
donnée de la simulation discrète, la classification des paramètres de contrôle et de leur domaine 
d'influence sur les comportements est également précieuse pour la compréhension des origines 
microscopiques des lois constitutives macroscopiques. 
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